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Probable detection of H 2 D + in the starless core Barnard 68 
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ABSTRACT 



Context. The presence of H 2 D + in dense cloud cores underlies ion-molecule reactions that strongly enhance the deuterium fractionation of 
many molecular species. 

Aims. We determine the HiD + abundance in one starless core, Barnard 68, that has a particularly well established physical, chemical, and 
dynamical structure. 

Methods. We observed the ortho-H 2 D + ground-state line lio-ln, the N 2 H + 7=4-3 line, and the H 13 CO + 4-3 line with the APEX telescope. 
Results. We report the probable detection of the o-H2D + line at an intensity r m b=0.22 ± 0.08 K and exclusively thermal line width, and find 
only upper limits to the N 2 H + 4-3 and H 13 CO + 4-3 intensities. 

Conclusions. Within the uncertainties in the chemical reaction rates and the collisional excitation rates, chemical model calculations and 
excitation simulations reproduce the observed intensities and that of o-H2D + in particular. 

Key words. ISM: abundances - ISM: individual objects: Barnard 68 - ISM: molecules - Submillimeter 



1. Introduction 

The densest and coldest cores of interstellar molecular clouds 
are receiving much attention as the possible precursors of 
star formation. Often referred to as starless or pre-protostellar 
cores, they are characterized chemically by a large depletion of 
molecules as they freeze out onto dust grains, and an associated 
increase in the relative abundance of deuterated isotopomers 
of numerous species (e.g.. iKuiper et al. 1996; Ceccarelli et al 
1998 [ lLoinard et alJl2002i Bacmann et alJl2003l ICaselli et al 
2003t lStark et al.l2004 . Phis increase is explained by the reac- 
tion H3 + HD H2D + + H2, which is favored in the forward 
direction at low temperatures, and subsequent ion-molecule re- 
actions involving H2D + . Detection of H2D + in dense cloud 
cores directly tests this proposed mechanism. Through its 
spectral line shape it also probes the core's velocity field, 
in regions where all othe r molecular tracers are strongly de- 
pleted. First predicted by Dalga rno et alJ ( 119731) . the ground- 
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state transition of ortho-H2D + , 1 10-1 11, lies in the submillime- 
ter wavelength range near 372 GHz, in a region with poor 
atmospheric transmission. Under excellent observing condi- 
tions on Mauna Kea (Hawai'i), HbD + has been successfully de- 
tected toward astronomical sou rces, including the young stel- 
lar object N GC 1333 IRAS 4A JStark et al.lll999l) . the starless 
core L1544 dCaselli et al.l 2003J), and the circumstellar disks 



of TW Hya and DM T au dCeccarelli et"al~ll2004l but also see 
iGuilloteau et al.ll2.QO6h . This Letter reports the probable de- 
tection of the H2D + 1 10-1 11 line toward the well-studied star- 
less core Barnard 68 conducted with the Atacama Pathfinder 
Experiment (APEX 1 ) telescope on Chajnantor (Chile). 

Barnard 68 (B68) is one of the most extensively studied 
starless cores. Using stellar extinction measurements, its den- 
sity structure has been foun d to be well matc hed by a near- 
critical Bonnor-Ebert sphere ( Alves et al. 2001). Bonnor-Ebert 
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1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration of the 
Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
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spheres dEberllll957t I Bonnorlll958l) describe the equilibrium 
configuration o f self-gravitat i ng clo ud cores just before the on- 
set of collapse. iBergin etaT, (2002) found strong depletion of 
C ls O in B68. lLada et all d2003l) analyzed the velocity field of 
B68, and found that the line widths in the center are close to 
thermal, leaving no room for significant turbulent motion. In 
the outskirts of the core, line centroid shifts suggest a non- 
radial pulsating motion. 

With the successful commissioning of APEX, the 372 GHz 
ground-state transition of ortho-H2D + has come within reach 
of regular observing, given the good local weather conditions. 
This Letter presents the probable detection of H2D + toward 
the center of this starless core obtained during the Science 
Verification of APEX (iJ2J and discusses the emission strength 
in the framework of a chemical model including depletion and 
deuteration (|3j. 

2. Observations and results 

The APEX telescope observed the H2D + lio-ln line at 
372.421385 GHz on 2005 July 24 and 25 using the APEX-2a 
receiver and the FFTS backend with a bandwidth of 1 GHz 
and 16384 channels. This frequency setting also covers the 
N 2 H + 7=4-3 line at 372.6725090 G Hz. The tele s cope w as 
pointed at the Ay peak measured by lAlves et all J200ll) at 
«2000 = 17 h 22m38 s .6 and <5 2 ooo = -23°49'46'.'0. The observa- 
tions were taken in position-switched mode, using an emission- 
free reference position. Pointing was checked on the nearby ob- 
ject RAFGL 1922. During the observations the source was at 
elevations of 25°-40°. A precipitable water vapor column of 
0.47 mm resulted in DSB system temperatures of 150-250 K. 
On 2005 July 25, the H 13 CO + 7=4-3 line at 346.998546 GHz 
was observed with a similar set up. After the observations, the 
velocity scale of the spectra was recalculated by the telescope 
staff to correct for a small (0.1 km s ) inaccuracy during data 
taking. 

After careful inspection of the data, the individual 30 sec 
scans were averaged using the CLASS software package; a 
total integration time (on+off) of 36 min was obtained for 
the H 2 D + line (12.6 min for H 13 CO + ). A sinusoidal baseline 
with a period of ~ 350 km s _1 was removed, followed by a 
first order baseline in the area surrounding the expected line. 
After smoothing the spectral resolution to 0.098 km s -1 , a rms 
noise level of 0.059 K on the T\ scale was found (0.058 K 
for H I3 CO + in 0.11 km s _1 channels), adopting a forward ef- 
ficiency of 0.97. The data were subsequently transformed to 
the main beam antenna temperature scale by division by 0.73, 
the recom mended mean beam efficiency in the 345 GHz band 
llGusten et all2006l) . 

The H2D + line was detected at a signal-to-noise ratio 
of ~3 (Fig. Fitting a single Gaussian line to the spec- 
trum yields r mb =0.222 + 0.082 K, V LSR =3.36 + 0.05 km s _1 , 
AV=0.33 + 0.1 km s _1 , and an integrated line intensity of 
0.078 + 0.015 K km s _1 . Although the line peak is detected 
at only 2.7<x, the integrated intensity is detected at 5.2<x, 
and we argue that the detection of H2D + in B68 is proba- 
ble (only deeper integration can make the result more secure). 
Fig.[3]in the on-line material shows the H2D + spectrum over a 



200 km s range. Out of the 1018 channels in this part of the 
spectrum, only 4 exceed the 3cr level (0.4%), as statistically 
expected. The significance of the detection is further supported 
by the close match to the source Vlsr of 3.36-3.37 km s _I and 
the purely t hermal line width of 0.33 km s _1 (H 2 D + at 10 K), 
as found bv lLada et 

ID feOOfl for N 2 H + and C ls O. All other 
peaks are much narrower (1-2 channels). Because H2D + re- 
tains a high abundance at the center of the core (see below), the 
apparent absence of turbulent motio n in the H?D + line provides 
strong support for the conclusion o flLada et al.l 12003 ) that B68 
is exclusively thermally supported. This situation is very differ- 
ent from the velocity field in, e.g., L1544, w hich shows signifi- 
cant v elocity gradients and wider H2D + lines Jvan der Tak et alJ 
2005). No detection was made of either the N2H + 4-3 line to 
a 2cr upper limit of 0.16 K or the H 13 CO + 4-3 line to 2<x of 
0.18 K. An emission peak at the correct Vlsr for the N2H + 
transition and the expected thermal width of 0.075 km s _1 is 
likely noise. 

In the documents accompanying the data release, it was 
noted that lines tuned in the upper sideband of APEX-2a could 
be too strong by 40%. This effect has been shown to be present 
for 12 CO 7=3-2, and is thought to be due to a resonance near 
the CO line frequency. It is therefore unlikely that our H2D + 
observations suffer from this effect, although no sideband-ratio 
measurements are available to confirm this. 

3. Discussion 

At the detected line strength of 0.22 + 0.08 K, the H 2 D + line 
is weaker than detections toward other starless cores such as 
L1544 (~ 1 K). Two aspects distinguish B68 from L1544. 
First, B68's H2 column density is lower th an that of LI 544, 
at 3.6 x 10 22 cm- 2 and (6-13) x 10 22 cm" 2 dAlves et al.ll200lt 
IWard-Thompson et all I2OO21 1 19991) . Secondly, B68's central 
density of 2 x 10 5 cirT 3 is lower by a factor of 6 than that 
of LI 544. Given that the critical density of the ground-state 
transition of ortho-H2D + is ~ 2 x 10 6 cm -3 (for our adopted 
collision rate; see below), a careful analysis is required to test 
whether the decreased line strength is due to (a combination 
of) sub-thermal excitation, lower total column density, or lower 
abundance of this key deuterated species. 

From the strength of the H2D + line, assuming optically 
thin conditions, we derive a beam-averaged column density of 
Af(H 2 D + )=(1.5 + 0.14) x 10 12 crrT 2 for an excitation tempera- 
ture of 10 K. This corresponds to thermal equilibrium excita- 
tion at the kinetic temperature in B68. Given the critical density 
quoted above, it is highly likely that the excitation is subther- 
mal, requiring a larger column density to reproduce the emis- 
sion. For example, an excitation temperature of 5 K implies 
N(H 2 D + )=(1.0 + 0.3) x 10 13 cm" 2 . This column density cor- 
responds to an average H2D + abundance with respect to H2 of 
(2.9 + 0.9) x 10- 10 . 

Further insight in the H2D + abundance can be obtained 
from modeling the chemistry in the B68 core followed by a 
statistical equilibrium calculatio n of the H?D + exc itation and 
line radiative transfer. Following Vastel et al. (2006) we use the 
simple chemical model of Casel li et alJ (|2002), updated to in- 
clude the multiply deuterated forms of HJ and the new val- 
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Fig. 1. (a) Detection spectrum of H2D + 1 k>-1 i i toward the cen- 
ter of B68 (histogram), with the best-fit Gaussian line profile 
superposed (thin line). The model of ^yields an identical pro- 
file for the modified collisional excitation rates as discussed, 
(b) Non-detection of N2H + 4-3 at the same position, (c) Non- 
detection of H I3 CO + 4-3 at the same position. 



ues of the binding energies of CO and N 2 llOberg et alJE ()05). 
B68 is modeled as a spherical cloud with radius 12,000 AU 
and the density profile f itted to the Bonnor-Ebert solution by 
iTafalla et all J2002i I20Q4I) . Th e temperature profile follows that 
of L1544 (lYoung et al.l2 004). resulting in a temperature gradi- 
ent from 13 K at the core edge to 9.8 K at its center. 

The chemical model starts with fully undepleted abun- 
dances of CO, N2, and O, and lets the species interact with 
and freeze-ou t on the dust gra ins which follow a MRN size 
distribution (Mathi s et al Jl 19771) with a minimum grain size of 
0.005 /im. Assuming that the ion chemistry is fast compared 
to depletion, the molecular ions follo w from the in stantaneous 
neutral abundances. Following Vaste feTall J2006I) a cosmic- 
ray ionization rate of 1.3 x 10~ 17 s" 1 an d sticking coefficient 
to S = 1.0 dBurke & Hollenbachlll983l) are adopted. Finally, 
for H2D + an ortho-to-para ratio of 1 .0 is adopted as expected 
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Fig. 2. Predicted abundances of several key species in the H2D + 
modeling (solid lines) and H2 number density (dashed line) as 
functions of radius in B68. 



at low temperatures fcrerlich et alJl2002h . We note that the re- 
sults of these calculations may depend quite sensitively on the 
chemical reaction rates and their dependence on temperature 
(Emprechtinger et al. in preparation), including the H2D + de- 
struction reaction with H2 dSchlemmer et alJl2005l and Asvany 
et al., in preparation). 

Figure |2 shows the resulting abundance profiles of some 
relevant species. Typical H2D + abundances with respect to H2 
at the center of B68 are predicted to be 4 x 10" 10 . The corre- 
sponding H2D + column density is 5.0 x 10 12 cm" 2 , close to the 
beam-averaged column density of 1 x 10 13 cm -2 found from 
the simple analysis above. 

With these abundance profiles, the statistical equi- 
libri um excitation is c al culate d using the code of 
iHoeerheiide & van der TalJ J2000I) . The velocity field is 
assumed to be static with only thermal broadening and no 
turbulent motion. We ignore the signi ficant motion in the outer 
regions of the core reported by Lad a et aD J2003h . since our 
observations are taken to the core's center and the observed 
spectrum shows no indication of line broadening. The resulting 
spectrum is convolved in a 17" FWHM beam appropriate for 
APEX at these wavelengths, yielding intensities on the T m \, 
scale; a distance of 100 pc is adopted for B68. 

These calculations require reliable collisional rate coeffi- 
cients, which are not a vailable in the literature. We follow 
van derTaketalJj2005l) in adopting estimated rates, and inves- 
tigate the effect of their inherent uncertainty. Using these rates, 
we find a H2D + 1 ^-1 1 1 intensity of 0. 10 K, lower by a factor 
of 2 than observed. Increasing the collisional rates by a factor 
of three, entirely within the estimated uncertainty, increases the 
emergent intensity to 0.24 K. Varying the rates by factors of 10 
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up or down produces a range of intensities of 0.011-0.48 K. 
The opacity of the line is 0.13 and the excitation temperature 
ranges from 3 to 6 K, consistent with the simple assumptions 
made above. 

We conclude that the observed intensity of H2D + is con- 
sistent with chemical model predictions, but that the uncer- 
tainty in the available collisional rate coefficients precludes 
a detailed comparison, especially for subthermal excitation. 
Determination of more reliable rates is clearly warranted. 
Ivan der Tak et al.l 12005 ) reach the same conclusion, including 
a preference for collision rates enhanced by factors 3-10; like 
these authors, we refrain from recommending adjusting the col- 
lision rates because of the uncertainties in the chemical model- 
ing. 

Using collision al rates from the Leiden data base 
( Sch oier et al.l 120051) . our model calculations are compatible 
with the upper limits on H 13 CO + and N 2 H + 4-3 r eported here, 
as we ll the N 2 H + and C 1 8 Q 1-0 detections of iBergin et alJ 
J2002I) and N 2 H + 3-2 from lCrapsi et alJfcOOSh . Quantitatively, 
we find 19 mK for H 13 CO + 4-3 (observed: < 0.15 K, 2cr); 
0.79 K km s 1 for C ls O 1-0 (0.85 K km s 1 ). Our model over- 
produces the N 2 H + emission by about a factor of 2; we find 
0.38 K for N 2 H + 4-3 (observed < 0.16 K, 2(f), and values 
of 3.8 K km s" 1 for N 2 H + 1-0 (observed, 2.5 K km s" 1 ) and 
0.35 K km s" 1 for 3-2 (observed, 0.17 K km s" 1 ). 

We conclude that our detection of H 2 D + toward B68 is 
consistent with models for the deuterium chemistry in starless 
cores which also reproduces other observed lines. However, 
the lack of reliable collision rates precludes any stronger state- 
ments about the H 2 D + chemistry. Given the pivotal role of 
H 2 D + in deuterium chemistry and the well studied nature of 
B68, calculation of such rates is warranted, especially now that 
observations of H 2 D + ground-state transition are possible with 
APEX. In the future, ALMA observations of the resolved emis- 
sion of H 2 D + and other species will provide further insight into 
the chemical state of B68, and only with reliable collision rates 
will quantitative analysis of these observations be possible. 
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Fig. 3. Detection spectrum of H2D + 1 io-l n toward the center of B68 over a baseline range of 200 km s 1 . Noise levels at + 
+2cr, and ±3cr are indicated. The arrow indicates the V L sr of B68 of 3.36 km s _1 . 



